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We explore changes in abundance patterns that occur late in the r process. As the neutrons 
available for capture begin to disappear, a quasiequilibrium funnel shifts material into the large 
peaks ai A — 130 and A = 195, and into the rare-earth "bump" a.t A = 160. A bit later, after the 
free-neutron abundance has dropped and beta-decay has begun to compete seriously with neutron 
capture, the peaks can widen. The degree of widening depends largely on neutron-capture rates 
near closed neutron shells and relatively close to stability. We identify particular nuclei the capture 
rates of which should be examined experimentally, perhaps at a radioactive beam facility. 



I. INTRODUCTION 

The r process, which synthesizes roughly half the el- 
ements with atomic number A > 70, proceeds through 
neutron capture and beta decay 0J|,p. Through most 
of the process, we know, capture is much faster than beta 
decay, so that very neutron-rich and unstable nuclei are 
temporarily created before disappearing as the process 
peters out. The question of where the nucleosynthesis oc- 
curs, however, is still unanswered. An evacuated bubble 
expanding behind a supernova shock waves is a promis- 
ing candidate for the site, but not yet the clear choice. 
The scenario is more convincing if the initial expansion 
is very rapid Q , but simple models of fast adiabatic ex- 
pansion 1^ imply that the neutron capture must finish 
in less than a second. In traditional simulations, where 
capture must often wait for a nucleus to beta decay, the 
process takes 2 or 3 times that long. 

Simulations from, e.g., ref. [|j demonstrate that the 
entire process can take place quickly if neutron capture 
populates nuclei farther from stability (and thus shorter 
lived) than usually thought. Despite initially forming at 
lower Z and A than traditional work suggests, the simu- 
lated abundance peaks end up at the right atomic num- 
bers. The apparent reason is that nuclei in the peak move 
up quickly in Z through an alternating sequence of beta- 
decay and neutron capture near the end of the r process, 
when the supply of free neutron begins to run out but be- 
fore neutron capture completely stops. But how do peaks 
maintain themselves during this late time, as beta decay 
drives each nucleus at a different rate towards stability? 
This question is actually more general than the rapid- 
expansion scenario; a quick move towards stability while 
the neutron abundance drops, though most dramatic if 
the path is initially very far away, in fact characterizes 
all bubble r-process simulations that produce something 
like the correct abundance distribution. And the question 
is linked to a broader issue: the significance of neutron 
capture once the process has slowed down so that it must 
compete with beta decay. What happens if the capture 
rates then are faster or slower than we think? In which 
nuclei do capture rates have the largest effects on final 



abundances, and can the rates there be measured? These 
are the kinds of issues we address here. 

The late stages of neutron capture and beta decay 
in the r process are much different from what precedes 
them. When the neutron-to-nucleus ratio R is much 
larger than 1, equilibrium between neutron capture and 
photodisintegration is a good approximation; the "path" 
consisting of the most abundant isotopes for each ele- 
ment Z is far from stability and moves relatively slowly. 
The term "steady state" is sometimes used to refer to 
this period, which ends when the neutron-to-seed ratio 
R falls below a few. To see what happens next, we note 
that the neutron separation energies along the path in 
(n, 7) <-> (7, n) equilibrium are related to the neutron 
number density n„ , which is proportional to R for slowly 
changing matter densities, by 
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where m„ is the neutron mass. When R drops below 
about 1, a nucleus on the path that beta-decays can no 
longer capture enough free neutrons to return to the path 
so that the path itself must move instead, inwards to 
higher neutron separation energy. The increased average 
neutron binding makes photodisintegration less effective, 
which in turn reduces the number of free neutrons still 
further, in accordance with eq. |l]. The process feeds on 
itself, causing R to drop exponentially and the r-process 
path to move quickly towards stability. Soon R becomes 
so small that (^,7) <-> (7,71) equilibrium begins to fail, 
and beta decay moves a good fraction of nuclei away 
from the path. Eventually, beta decay becomes faster 
than neutron capture and all remnants of (n, 7) <-> (7, n) 
equilibrium vanish. The inability of equilibrium to main- 
tain itself on the time scale of beta decay is usually called 
"freezeout" . 

To answer the questions about peak evolution and the 
significance of neutron capture at late times, we focus 
on two competing effects. The first, which dominates 
just as R falls below a few, when (n, 7) ^ (7, n) equi- 
librium still holds well, is a funneling of material into 
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moving peaks, most notably the small rare-earth peak 
lH, but also the larger peaks at ^ = 130 and 195. As 
time passes, the funnel fights an increasing tendency for 
the peaks to spread because beta decay and beta-delayed 
neutron emission compete harder with neutron capture. 
The interplay of funneling and spreading will imply that 
uncertain neutron capture rates, which are irrelevant as 
long as (n, 7) (7, n) equilibrium holds, become impor- 
tant fairly near stability, and should be determined there 
more precisely. 

We support our contentions with simulations of the 
neutron-capture part of the r process. In most, we ex- 
plore late times without worrying about a fast r process. 
We assume an exponential decay of temperature and den- 
sity with a relatively slow time scale of 2.8 seconds, the 
same as in ref. [Q. The distribution of seed nuclei with 
which we start is the post-alpha-process distribution of 
ref. [0, and we vary the initial temperature (and thus the 
initial density, since we keep the entropy equal to 300) 
and neutron-to-seed ratio. What we call our "standard 
simulation" is run under these conditions with the intial 
temperature fixed at Tg = 1.5. In section III, we worry 
more about a fast r process and therefore use different 
conditions. We do this in two ways: the first is to use 
use the same temperature and density dependence as de- 
scribed above, but a greater initial density (and thus a 
lower entropy). The faster neutron capture at high den- 
sity pushes the path farther from stability and leads to 
the formation of the A = 195 peak and the exhaustion 
of free neutrons in only tenths of a second. The second 
way is to use the conditions of Ref. where the tem- 
perature and density drop with a time scale of 50 ms and 
then level off at low values (Tg < 1, p < 100 g/cm^). The 
resulting drop in photodissociation rates again moves the 
path farther from stability, and the A = 195 peak forms 
in under a second. All our simulations use nuclear masses 
from ref. |^ and beta-decay rates from ref. |9|. 

The rest of this paper is organized as follows: In Sec- 
tion H, we discuss the action of funneling and spreading 
in the formation of the rare-earth element bump. Section 
m applies the same ideas to the large peaks, with empha- 
sis on the change in the peak's location and width as the 
path moves. The most important results appear in Sec- 
tion [V, which discusses neutron capture near the peaks 



and isolates particularly important rates that should be 
measured. Section ^ is a conclusion. 



II. FUNNELING AND SPREADING IN THE 
RARE EARTH REGION 

Ref. presented the basic dynamics of the near- 
freezeout funnel. It concluded that for much of the time 
when i? < 1, the system is still nearly in {n,j) <-> {"f,n) 
equilibrium, even as the path moves inwards. A kink 
at = 104, 106 soon develops (or grows stronger) when 
the path approaches a deformation maximum, which acts 



like a miniature closed shell. The nuclei near the bottom 
of the kink are further from stability than those at the 
top and thus have shorter beta-decay lifetimes (see fig. 3 
from Ref. which shows the path together with con- 
tours of constant beta-decay lifetime). The rate at which 
the path moves is governed by the average beta-decay 
lifetime along the path, which typically corresponds to 
a nucleus in the kink. The nuclei at the below the kink 
have shorter lifetimes than this overall average, and tend 
to beta-decay before the path moves, then capture neu- 
trons in an attempt to stay in equilibrium along the path. 
The nuclei above the kink have longer lifetimes and so do 
not usually beta-decay before the path moves, instead 
photodissociating to keep up with a path that is moving 
away from them. The neutron abundance is so low at 
this time that neutrons are essentially transferred from 
nuclei that photodisintegrate to those that capture. The 
net result is that nuclei both near the bottom and top of 
the kink funnel into it as the path moves, leading to a 
peak in the final abundances, whether or not one exists 
before R = 1. 

Another process, this one not discussed in Ref. acts 
to weaken the funnel. As R drops below 1, so does the 
rate at which neutrons are captured, since it is propor- 
tional to the neutron abundance. As a result, a nucleus 
in or near the kink will not always have time to cap- 
ture neutrons after it has beta-decayed; it may first un- 
dergo another beta decay and move away from the path 
of greatest abundances. Nuclei can emit neutrons follow- 
ing beta decay, moving them still further from the path. 
Thus, part of the growing peak begins to seep to lower 
neutron number N. Together with material from above 
the peak moving down in N, this process acts to wash 
out the peak in both N and A. 

Funneling and spreading counter one another, but as 
noted in the introduction, the two mechanisms reach 
their most effective points at the different times. Close 
to R = 1, when the path begins its inward trek, there 
are still enough neutrons so that spreading is slow and 
the funnel dominates. At very late times, by contrast, R 
is so small that (n, 7) <-!■ (7, n) equilibrium is seriously 
compromised and spreading is substantial. Eventually, 
neutron capture becomes slower than beta decay and the 
system freezes entirely out of equilibrium. After that 
capture essentialy stops and the only thing affecting the 
abundance distribution vs. A is some final spreading from 
delayed neutron emission. Beta decay without emission 
continues to move nuclei away from the equilibrium path, 
altering the distribution in N, but has no effect on abun- 
dances plotted vs. A (as they usually are). 

Our simulations make all these statements concrete. 
Fig. |l| compares the results of our "standard" simulation 
described in the introduction with the measured abun- 
dances in the rare-earth region as a function of A, show- 
ing the existence of a REE peak, in reality and in the sim- 
ulation. The simulated peak clearly relies on a kink that 
develops in the path because of the deformation max- 
imum but, as argued in ref. H, the mere existence of 
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a kink is not sufficient to fully produce such a peak; it 
achieves its full size only because of funneling. To see 
in more detail how the peak builds, we plot in fig. || the 
number of nuclei in three regions of iV — that just be- 
low the peak (iV=95 to 101), that including the peak 
(iV=102 to 106), and that just above the peak (A^=107 
to 113) — as a function of time for the run just discussed. 
The two vertical lines mark the points at which i? = 1, 
and at which beta-decay and capture rates are equal in 
the rare-earth region, causing the complete freezeout of 
(n, 7) <-> (7,n) equilibrium there. The bump develops, 
and then actually starts to shrink as material moves to 
higher A early in the process. But just before R = 1, 
as the path begins its inward move, it grows again. As 
noted above, photodisintegration is the dominant reac- 
tion above the bump and beta decay the important reac- 
tion below, so that material on both sides of the bump 
shifts inwards. After a tenth or two of a second, spread- 
ing begins in earnest; nuclei in the bump move slowly to 
lower N and material moves down from above to fill the 
trough above the peak, so that the abundance outside the 
peak starts to increase. To the right of the second verti- 
cal line, only beta decay (sometimes with the emission of 
neutrons) occurs, so that material is shifted downwards 
as the bump itself moves to lower N. 
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FIG. 1. Predictions of our standard simulation (solid line) 
for the final abundances of r-process nuclides versus atomic 
number A, and the measured abundances (crosses) scaled to 
the simulation. Note the peaks in the rare-earth element re- 
gion near A — 160. 

Although the brief region during which funneling dom- 
inates is evident in this figure, one may wonder whether 
the peak could form even if it was absent during the 
buildup phase along the path far from stability when 
R > 1. In ref. Q we argued that that was the case, and 
fig. H here provides more evidence. To make the figure 
we took the run discussed in figs. |^ and || and adjusted 
the abundances at i? = 1 so that the three regions of 
N were equally populated. We then let the run proceed 
starting from R ~ 1; the REE bump still formed. Fig. ^ 
clearly shows that funneling in the first tenth of a second 
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FIG. 2. Results of our standard simulation for the total 
abundances in the regions just below the peak (dotted line), 
in the peak (solid line) and just above the peak (dashed line) 
as a function of time. A drop in peak material is suddenly 
reversed when the neutron/nucleus ratio R nears 1, a point 
indicated by the first solid vertical line. The second solid 
vertical line indicates the time of freezeout in the rare-earth 
region. 



or so is responsible. [Some material is brought in from 
outside the range of the plot.] Here as before, later times 
show effects of spreading and the post-freezeout shifting 
of material downwards in N. 
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FIG. 3. The same as fig. except that the abundances 
have been adjusted when 7? = 1 so that the three regions 
all contain the same total. Now R = 1 occurs at f = and 
the solid vertical line indicates the time of freezeout. The 
formation of a peak here is due solely to funneling. 



III. FUNNELING AND SPREADING IN THE 
FORMATION OF THE A = 195 PEAK 

The fast r process of Ref. |^ relies on the formation of 
large peaks farther from stability than usually thought. 
Simulations show that the fast process works but but do 
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not explain how. Why should a peak that forms early 
at, e.g., N = 126, remain there when the path moves as 
neutrons are exhausted? In more traditional simulations, 
when the path is assumed not to move before freezeout, 
the usual explanation for peak buildup is approximate 
"steady beta flow" , which results in the longest live nu- 
clides building up the most. But this kind of buildup 
takes as least as long as the lifetime of the longest-lived 
nucleus, and the inward motion of the path we're dis- 
cussing here takes much less time. Something like steady 
flow therefore cannot be responsible for the existence of 
the peak at its final location in ref. ||]. What is? The 
answer is a tunneling phenomenon similar to that we've 
already discussed, though slightly more complicated be- 
cause instead of a kink we now have a long ladder of 
N = 126 isotopes populated at any given time. 

As already noted, the speed at which the path moves 
inward after i? < 1 is given by the average rate of beta 
decay along the path. This average rate tends to occur 
along nuclei in the middle of the ladders at = 82 and 
126. Thus the nuclei below them decay, and then capture 
into the peak as long as the funnel operates efhciently. 
A decay followed by a capture moves a nucleus up in 
Z, so material at the bottom of the ladder moves into 
the center of the ladder as the funnel proceeds. Further 
up the ladder, nuclei neither capture nor photodissoci- 
ate, since the path continues to run through these nuclei 
even as it moves toward stability. Instead they simply 
beta-decay, but more slowly than nuclei at the bottom of 
the ladder. The result is that the entire ladder shortens 
as the bottom moves up faster than the top. Above the 
ladder, for A^ just above 82 or 126, nuclei with slower 
beta-decay rates photodissociate into the peak, adding 
material just as in the REE region. These dynamics 
combine to move the peak up slowly in A at the same 
time as they heighten and narrow it. Though the large 
peaks clearly form during the steady-state phase of the 
r process, when i? >> 1, they are shaped and moved at 
later times as just described. 

These dynamics, however, can sometimes be masked 
by spreading. Whether or not they are depends on the 
temperature and density of the environment and on nu- 
clear properties. Fig. ^ shows the tunneling and spread- 
ing of material around the ^ = 195 peak in our stan- 
dard simulation as well as the two faster simulations de- 
scribed in the introduction. As in fig. |^, we plot sums 
of abundances in three regions of A^ — below the peak 
(A^ = 110 - 123), within the peak (TV = 124 - 126), and 
above the peak (A^ > 126) — as a function of time. In 
each case, the peak forms when i? > 1, but continues 
to grow for the few tenths of a second following R — 1. 
Much of this material comes from the photodissociation 
of material above the peak (dashed line in each plot). At 
later times, spreading takes over as material in the peak 
beta-decays to lower A'^. The onset and rate of spread- 
ing depends on how fast the neutrons are depleted after 
i? < 1, which in turn depends on the temperature and 
density. For a simulation at high density, as in part (3) 



of fig. H neutrons disappear very rapidly when i? < 1, 
and so in a short time beta decay and beta-delayed neu- 
tron emission win out completely over neutron capture. 
When the density and temperature drop quickly, as in 
part (2) of the figure, neutrons disappear more slowly. 
As a result, neutron capture competes with beta decay 
over a longer period of time, delaying the full onset of 
spreading. Nonevironmental factors can also affect the 
balance between tunneling and spreading. In anticipa- 
tion of our discussion of neutron capture in Section 4, we 
run simulations with the same three sets of conditions but 
a newer set of calculated neutron-capture rates, from ref. 
p2[ . Fig. H shows the tunneling and spreading of mate- 
rial in these simulations; the latter is less effective than 
in fig. ^. We expand on this point in the next section. 
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FIG. 4. Total abundances in the regions just below the 
peak (dotted line), in the peak (solid line), and just above 
the peak (dashed line) as a function of time for the three 
types of simulations discussed. Simulation (1) is the standard 
simulation, (2) uses the conditions of Ref. H, and (3) is the 
standard simulation at much lower entropy. 
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FIG. 5. Same as fig. but with simulations using a newer 
set of neutron-capture rates from Ref. | ]l2[ . 

The way the peak moves and gets shaped can be seen 
in fig. 0, which shows its time development in the three 
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types of simulations. While the peak initially narrows 
some as discussed above, it soon spreads, so that the ef- 
fect is barely visible. Fig. ^ shows the same development 
but with the newer capture rates. Here the narrowing 
of the peaks is evident, and is not erased by spreading 
at later times. The shifting of the peak to higher A is 
apparent in both sets of plots, and is most pronounced in 
the faster simulations where the peak forms much further 
from stability. 
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FIG. 6. The evolution of the 195 peak in late times for 
the three types of simulations, labeled as in fig. ^. The first 
frame shows the peak at i? 1, the third frame shows the 
final abundances, and the second frame is taken from a time in 
between, when R is less than 1 but much larger than its value 
at freezeout. For comparison, the dotted line in the second 
and third frames replots the abundances from the first frame. 
The scaled observed abundances are plotted as crosses. 
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FIG. 7. Same as fig. but with simulations using a newer 
set of neutron-capture rates from Ref. 



IV. SIGNIFICANCE OF CAPTURE RATES 

Funneling operates unhindered just a short time be- 
fore (n,^) <-> (7,71) equilibrium falters and spreading 
sets in. As we saw in the last section, once spreading is 
important, neutron-capture rates become so too. They 
determine how likely a nucleus that has beta-decayed is 
to return to the path before decaying again. Fast rates 
mean that (n, 7) <^ (7, n) equilibrium hangs on longer 
and spreading is delayed. Thus, the ultimate degree of 
widening a peak experiences depends on neutron-capture 
rates. To illustrate this point, we run simulations of the 
three types discussed above with four different sets of 
calculated rates |p|, ^ , mjl^ . These sets were calculated 
with different models for nuclear masses, slightly differ- 
ent treatments of the dominant statistical capture, and 
different assumptions about the importance of direct cap- 
ture. Not surprisingly, the rates can differ from one an- 
other significantly. Fig. ^ plots the ratio of the smallest 
to largest rates as a function of N and Z. When we 
use these rates in simulations (though all with the same 
mass model [||) we find variations in the final results for 
all values of A. We continue to focus on peaks, however, 
partly because the abundances are higher there than in 
neighboring regions, so differences are more significant, 
and partly because the differences in the left edge of the 
A = 195 peak are particularly noticeable. As we already 
saw in the last section, and as figs. ^ and |l^ show in 
more detail, the peak doesn't spread very much when 
rates are fast near the N = 126 closed shell. By contrast 
the slowest rates at these points cause the widest final 
peaks. These effects, incidentally, are particularly signif- 
icant for the Ref. [g| conditions, where (n, 7) <-> (7, n) 
equilibrium falters earlier because of the rapid drop in 
temperature and density, so that capture rates become 
important sooner. 
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FIG. 8. For each A'', Z the l og of the ratio of the highest 
neutron-capture rate in our set |2| Jl0yil| , [l2[ to the lowest. The 
darkest squares correspond to ratios greater than 1000, as 
indicated in the key. 

We can see the role of capture near the peak even more 
clearly by changing the rates only for between 123 and 
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FIG. 9. The four sets of neutron-capture rates, plotted as 
a function of for A'^ = 124 just below the closed neutron 
shell. 
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FIG. 10. Abundance curves from our standard simulation 
with the 4 sets of neutron-capture rates. Faster rates near the 
closed shell yield a narrower A = 195 peak (note the left edge 
of that peak). 

125. Fig. |l^ shows the results when those rates are multi- 
phed by 10 or 100, or divided by 100. When the rates in- 
crease, funneling becomes stronger and spreading weaker 
as (n,7) ^ (7,'^) equilibrium is partially restored. As a 
result, the final abundance peak at A = 195 narrows. 

It so happens that the nuclei at and just below closed 
shells are notoriously difficult to calculate ||l^. Com- 
monly used statistical methods may not be applicable 
for all those nuclei because of the low density of states 
at low energies [Q^. Rates of direct capture, which 
also plays a role, are uncertain because we don't know 
how much isovector dipole strength lies low in nuclei far 
from stability. To determine the astrophysical parame- 
ters in the r-process environment, it's therefore impor- 
tant to measure the rates in these nuclei where possible. 
Of course most of them are out of experimental reach 
for a long time to come. But the most important are 
actually relatively close to stability. 
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FIG. 11. Abundance curve when the capture rates for 
A'^ = 123 — 125 alone are changed. The solid line represents 
the results of with the rates from ref. the dotted line 
the results when those rates are increased by 10, the dashed 
line the results when the rates are increased by 100, and the 
dot-dashed line the results when the rates are shrunk by 100. 

To show why this is, we make even more selective 
changes, now for just 1 or 2 values of Z, in the rates 
for N between 123 and 125. Fig. |l| plots the root- 
mean-square difference between the abundance distribu- 
tion (within our standard simulation) when these rates 
are increased by 100 and and when they are unaltered, as 
a function of time. The ultimate degree of change depends 
strongly on which rates we change. Altering those below 
Z = 69 does little in the end because a) these nuclei 
are farther from stability, where the system is closer to 
(n,7) ^ (7,«) equilibrium and capture rates are nearly 
irrelevant, and b) Any changes that do occur have time 
to be diluted by spreading. Altering those above Z — 72 
doesn't do much to the final abundance pattern because 
the system has nearly frozen out of equilibrium, mak- 
ing neutron capture irrelevant because there are so few 
free neutrons and the temperature has dropped. The nu- 
clei for which changes do have large permanent effects lie 
between Z = 69 and 72 (Tm,Yb,Lu,and Hf), and corre- 
spond to the rough location of the path just before full 
freezeout, when neutron capture and beta decay com- 
pete on equal footing. Then dramatic differences in flow 
result from increasing the neutron-capture rates, differ- 
ences that are not erased by subsequent spreading. We 
see similar effects when the rates are decreased instead 
of increased. Fig. [l^ shows final abundance curve for a 
standard run with the rates of ref. [|| and standard runs 
in which the capture rates of the nuclei with Z = 69 — 72 
and TV = 123 - 125 are increased by 10 and 100. The dif- 
ferences are significant. All these statements remain true 
both when we make wide variations in the initial temper- 
ature, initial density, and time scale in our simulations, 
and when we use the conditions of ref. ||^. The reason 
is that in this range of Z, the calculated |^ beta-decay 
lifetimes of the N = 126 nuclei increase from about .07 
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s to 4.2 M. The path there slows down around the peak, 
giving the neutrons time to disappear through capture 
on nuclei in other regions. This happens whether in the 
time prior the path moves over a large distance (quickly 
at first, because of the fast beta decay rates far from 
stability) or over a shorter distance. Although it's theo- 
retically possible for freezeout to occur for Z < 69, the 
resulting peak would almost certainly be too low in A to 
match the observed abundances. 
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FIG. 12. Root-mean-square differences between the abun- 
dances with the rates of and with those same rates ev- 
erywhere except for a few nuclei with two values of Z and 
N = 123 — 125 (those are multiplied by 100), as a function 
of time. We used the standard-simulation conditions. Each 
curve corresponds to increasing a different set of rates. The 
nuclei with the largest effect on the final abundances have 
Z = 69 - 70 and 71 - 72. 

The rates for these few nuclei are thus the ones on 
which the final r-process abundances depend most sensi- 
tively, and measuring the associated cross sections would 
be useful; we could then better constrain the temper- 
ature and density during the r process. Unfortunately 
these nuclei are still far enough from stability that their 
cross sections may not be possible to measure, even par- 
tially through spectroscopic factors in transfer reactions 
with radioactive beams at RIA. A yield of about 10'*/sec 
is probably necessary for such experiments, while esti- 
mates of production at a RIA ISOL facility indicate 
that Z must be about 77 before yields will become that 
large. But we can approach the nuclei we're interested 
in, and see how measurements and calculations compare 
near the most critical region. 

We're more fortunate in the rare-earth region because 
neutron-capture rates are faster there than near the 
A = 195 peak, and freezeout therefore occurs closer to 



'^When we replace these lifetimes with the more accurate 
and faster ones calculated in ref. iQ, we find only slight dif- 
ferences in fig. nd. 
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FIG. 13. Final abundance curves in the simulation yielding 
fig. when only the rates of nuclei with A'' = 123 — 125 and 
Z = 69 — 72 are increased, by factors of 1 (solid line), 10 
(dotted line) and 100 (dashed line). Note the results of these 
simulations are nearly identical to those of fig. |ll|; the largest 
changes to the final abundance distribution are due to the 
modification of the capture rates of just these four nuclei. 



stability. Fig. |l^ shows what happens when we selectively 
change the capture rates of nuclei just below the kink, 
with N — 102 — 104, for particular values of Z. The nu- 
clei with the strongest effect now have Z = 62 and 63 (Sm 
and Eu). As before, the location of the most important 
nuclei is not very sensitive to initial r-process conditions. 
These nuclei are actually within RIA's reach. For the nu- 
cleus Z = 62,N ^ 102 (i64Sm), yields should be about 
lOVsec, for Z = 63, iV = 102 {^^^Eu) they should be 
about lOVsec, and for Z = 63, TV = 104 (^^'^Eu) about 
lO^'/sec |19]. Experiments to study their capture cross 
sections are worth considering. 
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FIG. 14. Same as fig. except for the rare-earth region, 
where rates of nuclei with A'^ — 102—104 and particular values 
of Z are changed. The nuclei with the largest effect on the 
final abundances have Z = 62 — 63. 

We have not discussed the A — 130 peak in any de- 
tail. Our conclusions there are more limited because we 
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do not reproduce the region below the peak very well. 
Abundances in that area are very sensitive to the out- 
com e of the alpha process, which we do not simulate. 



[12] T. Rauscher and F.-K. Thielemann, At. Data Nucl. Data 
Tables 75. 1 (2000); 



Ft p: ^ ^cjuasar . ph ysik . unibas. ch/ pub /tommy / astro /fits/ 



Non lcthclcss we do get a peak at A = 130 and gener- 



957fit_frdm.asc.g^ 



ally find that varying the N = 79 — 81 neutron-capture 
rates has the largest effect for Z = 48 - 51 (i27-i29cd, 
i28-i30in, i29-i3iSn, i30-i32Sb). RIA should be able to 
make enough of these isotopes to allow experiments. 



V. CONCLUSION 



[13] see S. Goriely, Phys. Lett. B436, 10 (1998), Astr. Astro- 

phys. 325, 414 (1997). 
[14] J. Engel, M. Bender, J. Dobaczewski, W. Nazarewicz, and 

R. Surman, Phys. Rev. C60 (1999) 14302. 
[15] See the RIA web page 



http://www.phy.anl.gov/ria/index.html 



For most of the r process, neutron-capture rates are 
irrelevant because they are fast enough to maintain equi- 
librium with photodisintegration. But at late times, the 
situation is different. Our investigation of funneling and 
spreading led us to identify particular nuclei relatively 
close to stability whose neutron-capture rates have sig- 
nificant effects on the shapes of peaks. It would be nice to 
have experimental information about these nuclei, even if 
indirect, e.g. spectoscopic factors through {d,p) neutron- 
transfer reactions. Measurements of the capture rates 
themselves in other nuclei closer to stability would also 
be useful; they would help tune models, which could then 
be better extrapolated to the important nuclei identified 
here. A full understanding of the r process would then 
be a little closer. 
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